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Abstract 

The total cross-section o~ tot and the jet cross-section o~j et differ at a proton- 
proton collision. The latter is divergent if arbitrarily small transverse mo- 
menta are allowed. Even with some fixed lower p± cutoff, o~j e t increases 
much faster than a tot at high energies. We have in this paper studied how 
the divergence could be tamed by colour screening effects among the par- 
tons. 

To do this we have built a proton model where we assign momenta, positions 
and colour charge to all partons in the proton. 

We find that the relative behaviour of the cross-sections can be better un- 
derstood by the inclusion of this effect. 



1 Introduction 

The Standard Model is the theory that all modern particle physics is based on. It tells 
you which particles that exist and how these particles interact with each other. The 
Model divides the particles into two main groups, matter particles and gauge bosons. As 
the name indicates, all matter is built up from the matter particles. The gauge bosons 
are the particles that mediate the interactions "forces" between the matter particles and 
themselves. The matter particles are divided in two different groups, the quarks and the 
leptons. These, in their turn, are placed in groups called families according to: 

(I) 




. . ( 2 ) 
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The upper ones are the quarks and they are named up, down, charm, strange, top and 
bottom. The quarks u, c, t all have (2/3)e in electric charge and d, s, b have — (l/3)e. 
All the matter that we have around us is made up of only the first family, that is up 
and down; the rest of the quarks exists today only inside accelerators and in some very 
energetic cosmic objects. 

The particles in Eq (H) are the leptons. The most famous of the leptons is the 
electron, but there are five more. The electron e, the muon \i and the tau r all have 
electric charge — e, while the neutrinos v are electrically neutral. 



Then there are the gauge bosons that mediate the different forces between the particles. 
There are four fundamental forces in the universe: 

Interaction Gauge bosons 

gravitation graviton 

electromagnetic 7 (3) 

weak W+ W", Z° 

strong gj, % — 1 . . . 8. 

Gravitation interacts with all matter, but it is so weak that it can be neglected almost 
always in particle physics. Then there is the electromagnetic force which is responsible 
for holding the atom together. The weak force is e.g. responsible for the /3-decay in the 
nuclei. In the standard model one has succeeded to combine the electromagnetic and the 
weak forces to one — the electro-weak force. And last there is the strong force that holds 
the quarks and nucleons together ]1J. 

In this paper we are going to look at the strong force and its gauge-bosons, the gluons, 
so a deeper study of this is in order. The strong charge is called colour-charge. Instead of 
+ and — there is red, blue, green and their anti-colours antired, antiblue and antigreen. 
All matter in the universe is built up of the quarks and the leptons. The quarks can't 
exist in isolation, but must be combined in colour-singlet states. This can be done in two 
main ways. You can combine a quark (colour) with an anti-quark (its anti-colour), this 
is called a meson, or you can take one quark of each colour (e.g. red-green-blue), this is 
called a baryon. The gauge bosons of the strong force, the gluons, can themselves carry 
colour charge. This means that the emission or absorption of a gluon can change the 
colour of the quarks as seen in Fig. [I]. The gluons have colour charge either as colour and 
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Figure 1: A gluon can change the colour of a quark 



a different anticolour (e.g.red-antiblue,r&) or are in colour neutral states (e.g.-^(rr — gg))- 

Almost all research in particle physics is done at huge accelerators. One is located 
near Chicago, called the Tevatron, and it is part of a research center called Fermilab. The 
accelerator ring is buried under the ground in a tunnel that has a circumference of about 
6 km. Here one collides protons (or actually protons and antiprotons) at energies as high 
as Ecu = 1-8 TeV. 

It is such collisions that we are going to study in this paper, or a bit more exactly, the 
cross-section of proton-proton collisions . The cross-section is a measure of how probable 
a reaction is. One can think of it as an effective radius of the incoming particles. If the 
two incoming particle's "radii" overlap each other we have an interaction. In our case we 
confront two different cross-sections, the proton-proton total cross-section Otot and the one 
called a j e t. In Oj e t we look at the particles inside the proton and their collisions instead. 
The particles inside the proton are usually called partons, and these include both quarks 
and gluons. It is these collisions between the partons that one can calculate theoretically. 
From these calculations one can create a theoretical a tot that you can compare with the 
experimental total cross-section. The measured cross-section increases much slower with 
the energy than the theoretical. There already exists some explanations to this, such as 



multiple interactions, see Section |272| , but these are not enough to explain the gap between 
the experimental data and the theoretical predictions. 

What we have done is to build a simple model of a proton with all its particles, 
partons. We have for all these partons calculated the momenta, positions and the 
colour-charges. We have also constructed the model so that the whole proton must be 
colour-neutral. What we are hoping to see is that the partons will screen each other and 
thereby the cross-section will be reduced. If the energy is high, the number of partons is 
large and this screening effect should be bigger. That is, the increase of the cross-section 
should be slower, which is exactly what is desired. We show that such a screening indeed 
appears to be an important physical mechanism. 



This paper is organized as follows: We will in Chapter 2 make some theoretical ar- 
guments about the evolution of our proton model and discuss the screening effect. In 
Chapter 3 we will describe our model for the proton. Later in Chapter 4 we uses this 
model to observe how a gluon interacts with a proton. And last in Chapter 5 we summa- 
rize and give some conclusions. 
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Figure 2: A proton-proton interaction where two partons scatter against each other 

2 Theory 

If we collide two protons with each other, see Fig. || , we have a parton-parton cross-section 
according to: 

da jet = J dxJixuQ 2 ) J dx 2 f(x 2 ,Q 2 ) J ^ d P±- ( 4 ) 

The first integral factor in Eq (|4]) expresses the flux of partons from proton 1. The 
second term is then the flux of partons from proton 2. The last term is a measure of the 
actual probability that two partons will interact with each other. 

The parton function f(x, Q 2 ) expresses the probability to find a parton with momen- 
tum fraction x of the incident proton momentum, while Q sets the virtuality scale at 
which the proton is probed. Q is often in hadron-hadron collisions represented by p± 
since these variables scale approximately the same, that is 

do-j e4 = J dx l f{x u p\) J dx 2 f(x 2 ,p 2 ± ) J -^T d P±- ( 5 ) 

The parton functions f{x,i ) p 2 A _) in Eq (|5|) increase very fast at small x-values (see Sec- 
tion 2~I for further details), that is small momentum for the partons. They also increase 



for bigger virtuality, Q. This will eventually lead to a parton-parton cross-section that is 
bigger than the whole proton-proton cross-section. 
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2.1 Evolution Equations 



When you are studying how a hadron, e.g a proton, is built up you have to do this by 
sending a particle, a probe, at the hadron. This probe then scatters against a parton 
inside the hadron and by this you can draw conclusions about how the hadron is made 
up. If the probe has a large Q (is very virtual) it can resolve more details and we see 
more partons. 

That there are details to resolve is a consequence of the quantum mechanical fluctua- 
tions. Virtual fluctuations a —>■ be —>■ a appear and disappear continuously, with lifetimes 
related to the virtuality scale according to the Heisenberg uncertainty principles. The 
fluctuations can be nested, thereby giving rise to whole cascades of fluctuations. 

There are three different processes that can happen in the parton cascade. These are 
the three vertices of the strong interaction of the first order. A quark can split into a 
quark and a gluon. And a gluon can either split into a quark-antiquark pair or into two 
gluons, see Fig. [| These processes are of course reversible, so e.g a quark can absorb a 
gluon. 




Figure 3: The 3 main vertices in the strong interaction 

The daughters share the momentum of the mother between themselves according to 
some probability distributions. It is e.g. a large probability that the gluon emitted by 
the quark has very small momentum. These probability distributions are described by 
the Altarelli-Parisi splitting kernels [pi 

4 1 + z 2 

*W*) = 3^7, (6) 

p„ w . 3<*!)! m 



where rif is the number of allowed flavors, usually 3-5. The variable z is the fraction 
of the mother momentum that the first daughter takes (i.e. the quark in Eq @). The 
other daughter then gets (1 — z) of the momentum. Both in Eq (|6]) and Eq ([?[) there are 
singularities, in Eq @ at z = 1 and in Eq (0) both at z — 1 and z = 0. The physical 
explanation is that there, at any given time, is a cloud of gluons with very low momenta 
around a quark. In reality even if there are are infinitely many gluons so are their total 
momentum finite. One usually introduce a cutoff, z m i n , and branchings that includes 
gluons below this z m i n (or over (1 — z min )) are then not allowed. 
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The probability that we will have a branching a — > be at all is given by the evolution 
equation, the DGLAP equation ||: 

fd,2^W«) (»> 



d(hiQ 2 ) y 2tt 
or, if one introduces the evolution parameter t = In ^2 



dPa-bc f , a 8 (t) D . v , in , 
- dz— — P a _ bc (-2)- (10) 



dt y 2vr 

That is, for a small increase in t, dt, there is a probability dP that a branching a — > be 
will take place. So the evolution equation evolves in the virtuality, Q. This does not 
mean that the individual partons evolve in Q, but it is simply an easy way to pick the 
right value for the branchings at which the partons are resolved. 

Another way of expressing this is an inclusive picture; you could calculate the prob- 
ability of finding a parton at a certain x value. This becomes more like a mean-value of 
the distribution, but it is a picture that is often used since it is these distributions that 
you can measure at experiments, and not how a single parton develops. We have, for the 
inclusive picture, the following expression: 

dfi(x) _ f 1 f 1 a 



d(lnQ 2 



/ dyf)(y) f dz£Pj^(z) 6(yz-x). (11) 
Jx Jo Ztc 



The first integral in Eq (|TT|) simply says that the mother, j, has to have at least x 
momentum. The other integral one recognizes from our earlier expression in Eq (H), and 
the delta function ensures that the daughter, i, obtains x in momentum. If we use the 
delta function to eliminate one of the integrals we get 

mx) - rldy -fM^J-). (12) 



d(hiQ 2 ) Jx y JJX ^2n 

So this expresses number of daughters that will get x momentum from a mother of mo- 
mentum higher than x. But these partons with x momentum can also decay, so we have 
also an out going flow from this point that has to be subtracted from the expression. 

Once you have a parton distribution at some Q-scale you can by these evolution 
equations iterate towards higher Q and get the whole distribution above this Q. This has 
been done by many groups [ffl |5| for different starting values. One feature that is common 
for all groups is that the distributions are very high at small x-values. 



2.2 Multiple Interaction 

If two protons collide with each other it can happen that more than one pair of par- 
tons collides inside the protons. This we call multiple interaction 0. There are some 
experimental data that indicates this phenomena |], §. 

The total cross-section, a tot , is of the order 40 — 70 mb see Fig. ^. It is often parame- 
terized according to [|]: 

a tot = As"** + Bs e , (13) 
where s is the E C m squared, 77 ~ 0.5 and e ~ 0.08. 
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Let us, on the other hand, look at the parton-parton cross- sect ion, and make some 
rough estimates such as: fi(x,t) ~ x ( >~ 1 ~ e >, and introduce an x min to avoid to get singu- 
larities at x = and finally a p± m in- P±min is introduced as an integration limit that you 
after calculating the integral should put to zero. If you do all this you will get something 
like: 

CTjet = J &Xif{Xi,Q 2 ) J dx 2 f(x 2 ,Q 2 ) J ^-^P± 

(14) 



2 1 e 
1 S 



Pi nun PI nun 



As is clearly seen, one must have p± m in 7^ 0. This contradiction can, we hope, be explained 
by the screening effect, see the next section. If you do these calculations a bit more 



carefully, or uses a simulation program such as PYTHIA 6.1 [10| you get curve according 



to Fig. |]. One can increase p± m in but this will only move the problem to another region. 
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Figure 4: The total cross-section a tot and the parton-parton cross section <Tj et . P± m in is 
here set equal to 2 GeV 

The huge difference between the two curves can, in part, be explained by multiple 
interaction since we have: 



O'tot = 



a 



jet 



&0jetpair ^ljetpair &2jetpairs ' ' ' j 
&ljetpair '^•^2jetpairs 3(J3j e tpairs 



(15) 
(16) 



where we have to multiply by the number of jets in <jj e t because (Jij etp air tells about the 
probability to have i jet pairs. If we have i jet pairs we always has i such terms. This leads 
to that at an event with many jets (high energy) the (Tj et will be higher than <j tot . 
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This seems to explain the curves well, but unfortunately it's not enough. We believe 
that the necessary last correction can be explained by the screening effect, that is that 
the quarks will not interact as if they are free particles which has been assumed in the 
perturbation calculations of the aj et cross-section. 

2.3 Screening effect 

All hadrons are in colour singlet states. This means that if you can't resolve a hadron it 
will not interact strongly. That is if you e.g. send a gluon with a very long wavelength 
at a proton, the gluon will not be able to see the partons inside the proton, but only the 
proton as a single object. Since this object is in a colour singlet state the gluon will have 
problem to interact with it. On the other hand, if the gluon is very energetic and has a 
short wavelength, it will resolve much more partons and the screening effect will become 
much smaller. 

This can be studied by looking at the ratio A between the incoherent and the coherent 
sum of colour charges in the cascade. 



A = 



\Ek=iQke ipXk \ 


2 


2^k=i 


Qk\ 


2 



(17) 



Here is n the number of partons in the proton, qu is the colour charge of the fc:th 
parton and is the position of the fc:th parton. The gluon has a momentum of p. If the 
gluon has a low p the wavelength is so long that it won't matter where in the proton the 
partons are, see Fig §, and all terms will cancel. Or you could see it like this: if p — ► 
then A — ► (since Qk — 0) and if p — > oo then A — > 1 




Figure 5: The screening effect; A gluon interacts with the partons in a proton and de- 
pending on the wavelength it interacts differently. 
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3 The Proton-Model 



We start, in our model of the proton, at a very low virtuality scale Q. At this low Q we 
assume that we have a finite number of partons, or more specificly, that we have 3 quarks 



and 2 gluons [llj]. These 5 partons are then evolved to higher Q by Eq (|T0|) to some 



Qmax- As we do this we continuously assign each new parton with momentum, position 
and colour charge according to some given probability. We have decided the transversal 
position of the partons while in the momentum space the longitudinal momentum is 
selected. We do this since it allows us to build a model and still be consistent with the 
Heisenberg Uncertainty Principle. To have a colour-neutral proton we, to begin with, 
have a colour-neutral start configuration, but then as we evolve we keep all partons. To 
avoid getting too many partons we had to introduce a x m i n and branchings in which a 
parton had a momentum under this were not allowed. This rather abrupt cutoff did not 
change the distribution significantly, as we shall see. 



3.1 Start configuration in momentum space 

We choose to start with 3 quarks and 2 gluons. In this choice we were guided by other 
works done earlier. In the original in 1976 [0 one started with only the three valence 
quarks, and this was not enough to fit the data with experiments. The GRV group 
JTT| , |T3"| later made the assumption that you had to introduce some gluons and seaquarks 
already from the beginning, though only very few, and this seems to be sensible since 
they got a good fit to the experiments . We cannot use the GRV distribution itself, 
since they use a inclusive picture, while we need an exclusive picture to be able to assign 
momenta, colour charge etc. to all partons. 



In our model, each parton was given a fraction x of the total momentum from a 
distribution that was fitted to give values similar to what the GRV group had. What we 
are trying to get is a qualitative picture of the parton distribution, not to compete with 
GRV for the best configuration. We made the following ansatz: 

f(x) = Nx a {l-x) (3 , (18) 

where N is a normalization factor to get the right number of partons and a and (3 are 
the free parameters that we used to adjust the distribution. 

To get the proper normalization X^=i x« = 1 we made the following operation: 

I \ — — 

\X i ) norm ^5 

Z^j=l Xj 

This means that the original distribution is pushed towards the middle 
had to choose a < 0. Some trial and error finally gave that: 

a = -0.4 
/3 = 1.2 

Though even this did not give a very good distribution. We had too much partons at 
high momentum and too few at low, see Fig. || This we solved by starting our evolution 
at Q = 0.44 instead of Q = 0.48 as they had done. That means that when we were at 
Q = 0.48 some of the partons with high momentum had split to two other partons with 
lower momentum and we got a good fit, see Fig. |. 



(19) 

As it were we 
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Figure 6: Different start configurations for the momentum. Theoretical curve from [11 



3.2 Start configuration in coordinate space 

Our model is 2- dimensional. This is not a strong limitation, though, since the protons 
have a very high speed and thereby are strongly Lorentz-contracted. That is, they can 
with good approximation be taken to be 2-dimensional objects. And since the only thing 
that we are going to do is to send plane waves in the transverse direction of the beam at 
the proton, see Fig. ^ we can project the points onto the x-axis. So we put x = rcos9 
and chose 9 at random for each new branching. 



y 




Incoming plane wave 



Figure 7: The projection of the position 

We have introduced a characteristic radius, Tq, of the proton at r$ = 0.7 fm. The 5 
original partons are then randomly chosen according to some distribution. We tried some 
different distributions: 

• Gauss : /(r)dr oc e 2r ° rdr 
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• Exponential I: /(r)dr oc e r ordr 

r 

• Exponential II : /(r)dr oc ^—^-rdr 

Here is ro the characteristic radius. The difference in the final results of the parton shower 
did not depend much on which of these distributions we choose. We decided to take the 
Gaussian distribution. 



3.3 The evolution of the momentum 

We have chosen to work with an exclusive picture, that is we follow the development of 
each parton in the proton. This picture is easy to work with since you can iterate from 
a start configuration of partons. You know the probability for a parton to split up at a 
little step, dt, of the evolution parameter. So by using random numbers to simulate the 
probability of the splitting you can reproduce the actual process. In the exclusive picture 
we had: 

dPa^bc /"i-w a s {t) 

—7— = / dz— — P a _bcU), (20) 

with 

127F - 1271 

as ~ (33-2n / )ln(Q 2 /A 2 ) ~ 25? ( ' 

where n/ = 4 and t = ln(Q 2 /A 2 ). Since we have that z min = x m i n /x we can't have 
any partons under x m i n . For the quark channel q — > qg this gives, with help of the 
Altarelli-Parisi splitting kernels to decide -P q ^ qg , that the evolution equation looks as: 



dt 25t J z m i n 3 1 — 2? t 

You can, by simply putting in the different splitting kernels, get the evolution equations 
for g — > qq and g — ► (7(7. For a gluon you have to decide which of the channels it will take: 
either the gg or the qq. This is randomly chosen according to the different probabilities 
of the channels. 

The number of partons, N, at a given Q is then given by 

One can draw a parallel here to the exponential decay of a radioactive nucleon. It has 
a certain probability, P, to decay, but if you look at a large number of nucleons the 
probability that you will have a decay in a little time interval decreases as the time 
elapses since you have fewer nucleons left which can decay. By solving Eq ([23|) you get 
how many partons, N, you have at a given Q-scale. 

In this way you evolve your parton distribution until you reach Q = Q max . We have 
also introduced a x m i n to keep the number of partons in our model from exploding. All 
branchings in which any of the daughters has a momentum less than x m i n are forbidden. 
We have compared our model with values from the GRV group |IT], see Fig. |8] and 



Fig. H We have an excellent agreement if we include all partons, as seen in the upper 
curves . The antiquarks though are somewhat suppressed in our model, especially at 
low Q. This is not so strange, since we in our model assume that we don't have any 
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antiquarks in our start configuration and GRV have a mixture of them. At higher Q 
these differences shrink and the two curves approaches each other, see Fig. |^ . We have a 
fairly good agreement at low x-values also, see Fig. [II], and this is a bit surprising since 
we have introduced cut. 
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GRV, all partons at (r=1 
model, all partons at Q =1 
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Figure 8: The evolution of the momentum at Q 2 = 1. At Q 2 = 1 q is too low in our 
model. The quark distribution can be attained by subtracting the gluon and the antiquark 
distributions from the total. 
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GRV, all partonsat (r=10 
model, all partons at Q =10 
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model, gluons at Q =10 
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Figure 9: The evolution of the momentum at Q 2 = 10. The q curve of our model 
approached the theoretical curve. 
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Figure 10: The evolution of the momentum at Q 2 = 10 in log scale. x m i n = 10 
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3.4 The evolution of the position 



The distance a parton will reach from its origin of creation depends on a number of things. 
In our model we assign each new pair of partons with a distance, x, according to: 

R, 

x = Xq ± — cos 6 a/1 — Zi (24) 

Here is Xq the position of the vertex where the two partons were created and R is a random 
number between and 1. 

The first factor ^ comes from the Heisenberg Uncertainty Principle. A parton with 
high virtuality, Q, lives shorter than a parton that is almost real. The random number, R, 
is there since a parton has freedom to decay at any point between xq and forward to our 
maximum distance 4. We have experimented with this factor and chosen some different 
distributions. One, perhaps more physical, is that instead of choosing a random number 
linearly to chose them from a exponential distribution. This agrees with a quantum 
mechanical picture, it is a high probability that it decays early but there exists a chance 
that it will live very long. 

The second factor in Eq (p31) cos 8 is simply the projection of the partons position to 



the x-axis as discussed in section 3.2. The angle 9 is randomly chosen. 



The third and last factor comes from the kinematics of the system, according to Fig. |11 
We have, from the conservation of energy and momentum, the following expressions: 

{E + p g ) 1 = z(E+p g ) (25) 
(E + p z )(E - Pz ) = m 2 + p 2 ± (26) 

Here p z is the direction of the incoming particle and p± is the direction orthogonal to this. 
If we apply Eq (p6|) on the initial state (0) and the final state (1 and 2), we get 



= -Q 2 +pI v\ = -Q 2 +pI v\ (27) 

(E+p^ (E+p z ) 2 z(E + Pz ) (l-z)(E + Pz ) - { ] 



In the last equality Eq (^5j) was used. From this expression it is straight forward to get 
a relation between p± and z according to: 



p± = VT^ Q. (28) 

If a parton gets a large share, z, of the momentum it will not be able to deviate much from 
the original partons direction. Therefore the distance it will travel according to Eq (|2"1| ) 
is small. We have also studied results without this factor, and find that it is not critical 
for the qualitative results. 




Figure 11: The kinematics of a vertex 
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3.5 The colour charge 



Since the colour charge has 3 different charges (colours), instead of as in electromagnetism 
2, you have to introduce one more dimension to represent the colour space. That is, you 
will get a colour-plane, Fig [12], where the different charges of the particles are added 
to each other as vectors. Though this is not completely true, since we have the same 
phenomena here as in the spin-space. Even though the z component, S z , of the spin is 
zero the total spin, S, can be a finite value, since S z is the projection of the spin vector 
on the z-axis. The same thing happens in the colour space, we can have a colour neutral 
particle that is not a colour singlet. An example of this are the two colour neutral gluons 
(e.g.-^(rr — gg))- We don't make this distinction in our model though, and we treat the 
colour charge as if it simply were vectors in a two dimensional space. So, e.g. you get the 
gluon gb charge by adding the vectors of the green charge and the antiblue charge, and 
so on. 




Figure 12: The representation of colour charge in a 2-dimensional room 

This representation of the colour charge is not perfect, we have only 6 gluons, as is 
seen in Fig. O. The colour neutral gluons (2) are not represented in this picture, since 
they would then coincide with origo which would give a set of none-interacting gluons. 
Another thing that is not entirely correct is the length, or charge, that the gluons get. In 
this representation the length of the gluons is while it should be = \J~ifz = f ■ 

For our simple model of the proton these small errors are not very important and this 
representation is adequate. 

We generate the start configuration of the 5 partons as if we originally had only 3 
quarks that were colour neutral, e.g. one quark of each colour. These emits 2 gluons 
randomly chosen (though not gluons that are colour neutral), according to q — > qg. 

A possible start configuration could then be that the red quark emits a gluon with rb 
and the blue emits a gluon with bg~. And we have 1 blue quark, 2 green quarks, 1 gluon 
with rb and another with bg. If you add these together they will exactly cancel each other 
and we have a colour neutral configuration. 
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As we evolve our parton shower we follow the rules of QCD in each vertex. Whenever 
there is freedom to choose colour is this done randomly. As long as we follow these rules 
we are guaranteed to have a colour neutral state. However, if we keep all partons, and 
evolve all this we will get several hundreds of partons in our model and this is not desired. 
To avoid this we have introduced a x m i n , as mentioned before. Branchings which produces 
partons with x < x min are not allowed, but we keep all the partons in the parton shower 
to keep the colour charge to be zero. 

3.6 Number of partons 

Here we present a graph showing how the number of partons increase with the virtuality, 
Q, for two different energy scales. We adjust the energy by changing the variable x min 
and the lower x m i n the higher the energy. 
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Figure 13: How the number of partons increases with the virtuality, Q. 



3.7 A summary of the proton model 

Our Model of the proton is made in an exclusive picture. We start at a low Q with 
only a few partons and let these partons evolve in an iterative process to some Q ma x- To 
reduce the number of produced partons we had to introduce a cutoff in the momentum 
space. The branchings that produces partons with a momentum under this cutoff were 
not allowed. Despite this rather abrupt cut in the evolution of the parton shower we got 
a good distribution of partons and their momenta. We have compared our results with 
GRV in Fig. [8] to Fig. [l(] and the total parton shower coincides almost exactly. 

We will now use this model to study some different phenomena, but mainly the screen- 
ing effect. 
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4 A study of the Model 



The primary goal of this section is to study how our model of the proton interacts with 
a gluon, but also how the radius of the proton varies with the energy of the system, 



section |4.5| . When we are sending a gluon against our proton will we see that the screening 
effect comes into play and decreases the interaction. We are going to try some different 
start configurations and see how large the screening effect is, and draw some conclusions 



of this, section 4.4 



4.1 Introducing the Interference equation 



As we discussed in sec |2.3| is the interesting equation when we are studying the interaction 
of an incoming gluon against the proton the following: 



.4 



EjLi Qke 



2^fc=i \Qk\ 

Since we in our model have assigned both position and colour charge to all partons 
we can readily use this formula and study how A varies when the incoming gluon has 
different transverse momentum. According to the theoretical predictions we should have 
a curve that at low momentum for the gluon is zero and then gradually increase to the 
maximum value 1 at very high momentum. This is exactly what we have, see Fig. 14. 
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Figure 14: The resolution of the proton increases with the gluon momentum and therefore 
we approach an incoherent state. x m i n = 10 -3 
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4.2 A running Q 



In chapter |4.1| we did the following: We first evolved the parton cascade according to the 
evolution equations. Then we sent gluons with increasing momenta at the now evolved 
proton. 

Another, more physical, approach would be to probe the proton with a gluon that 
had approximately the same value on p as the partons has Q. As we mentioned earlier 
(chapter |2|) , is Q ~ p±. So this procedure is built upon the principle that a probe 
can resolve best at its own wavelength (or momentum). In practice this is done by 
continuously, as we evolve our parton cascade to higher Q, sending in a gluon that has 
the same p as we have Q at this point. 



0.8 



0.6 



0.4 



0.2 







Q follows p 
Q follows p above p=1 




0.5 1 1.5 2 2.5 3 3.5 4 

P gluon (GeV) 



Figure 15: The Interference function A(p) when Q follows the p scale. In the dashed 
curve is parton shower frozen at Q — 1 and downward to avoid the peak at p = 0.44. 



10" 



The peak at p = 0.44 in Fig. [15] on the upper curve should not be there. It arises from 
the fact that when we start our evolution at Q = 0.44 we have only 5 partons and these 
5 partons have no tendency at all to cluster, unlike the partons produced later in the 
evolution. Since the physics in this low region is little understood our start configurations 
can not be checked correctly and it is not certain that the ansatz with only few partons 
is applicable. What we have done is to take our distribution at a higher Q, in our case 
at Q = 1, and use this distribution all the way down to Q — 0. This is similar to what 
is often done to calculate jet cross-sections with parton distributions not defined below 
some Q 2 . If we do this we get the dashed curve in Fig. and the peak has disappeared. 
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4.3 Introducing a variable x m i n 

We introduced a cut in the momentum, x m i n , and partons under this limit was frozen. 
We did not discuss how this cut was chosen, nor its physical meaning. 

If we look at two incoming partons with a momentum of X\ resp. X2 we have: 

s = X1X2S, (30) 

where s is Eq M and s is the energy of the subprocess squared. p± can at the most be half 
of the energy of the subprocess (since we have two partons) and this will mean that we 
get: 



P± < 'I = xix 2 ^ p\ ~ x 2 min ^, (31) 



that is, we have 



x min — • (32) 

This is a somewhat crude calculation, but we see two things here. First, and not so 
surprising, that the higher total energy we have the better we can resolve partons with 
small momentum, i.e we get a lower x min . The other thing is that, for a fixed E C m, x min 
scales linearly with p±. So, as we increase the virtuality scale, Q, we should continuously 
increase . 



< 
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Figure 16: Introducing a variable x m i n . In the dashed curve is Q running with p while in 
the other curve is Q = 2. We have also included the curve for a fix x m i n = 10~ 3 . 

In Fig. [16] we have that x m i n = lO -3 ^^. This factor 10~ 3 can be varied and the curves 
will be suppressed for smaller factors. In the dashed curve with a running Q have we, as 
said in chapter fO| , frozen the parton distribution at Q — 1 and downward. As can be 
seen in Fig. [16] this effect is small. 
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4.4 The Screening Effect 

All the graphs in chapter 4 are a measure of the screening effect since they all tell you 
how well the gluon can resolve the partons inside the proton. By parameterizing these 
curves by some function S(p±) you would see the strength of the screening effect and how 
it increases with the energy. As a reasonable guess, we put: 

S (p±) = 2 P } 2 . ( 33 ) 
P± + P±o 

where p±o is a free parameter. As long as p± is large is S equal to 1, it is only when p± 
approaches zero that our correction term comes into play. 
In the expression for the cross-section <jj et we had 

do-j e t = J dx x f{x u p\) J dx 2 f(x 2 ,p 2 ± ) J ^r d P±j (34) 
or, if we rewrite it a little 

dx 1 f(x 1 ,p 2 L ) I dx 2 f(x 2 ,p 2 ± ). (35) 



dcTjet _ da 

2 ~~ A„2 



dp]_ dp 

It is the first term in this expression that has a singularity as p± — > since it goes 
like l/pj_. By introducing S(p±) (squared since we have two partons) into Eq fl3"5| ) this 
singularity disappears and we get a finite value: 

dcr^ „, sQ da 



S{p±) 2 -^ J dx 1 f{x 1 ,p 2 L ) J dx 2 f{x 2 ,p 2 ± ) 

dx 1 f(xi,p 2 L ) J dx 2 f(x 2 ,p 2 ± ) 



dp\ dp 

p\ da 

(P± + P±o) 2 d P± 

~ (p2 +p 2 y J teif&hPj.) J dx 2 f(x 2 ,p 2 ± ) (36) 

When introducing the screening function S(p±) we adjust the parameter p± e.g., so that 
S(p±) coincides with our curve at A = 0.5, see Fig. [H]. By doing this for different energies 
we can plot how the screening effect varies by the energy. We change the energy of the 
parton shower by changing x min since we had the relation that x min ~ 1/Eqm- 

The curve that has a fix Q can better be fitted with our screening function than the 
one with a running Q, see Fig. |I7| This indicates a shortcoming either of the model or in 
the ansatz of S. We will use a fix Q mostly, to allow a more sensible comparison. 

By plotting p± versus 1 / x min we get a good measure of how the screening effect varies 



by the energy. We have, in Fig [18] done this for some different start configurations and 
also some different evolution conditions. Since we do not know how much the partons 
separate themselves in a branching we will also look at events where the splitting is twice 
the size of Eq ([24]) 

1. A fix Q = 2 and a branching according to Eq ( F24"D 

2. A fix Q = 2 and a branching twice the size of Eq (p4j) 



3. A running Q and a branching twice the size of Eq 

We see in Fig. [18| that pj_ increases as the energy of the system is increased. By this 
we get that aj et increase slower than for a fixed p± and we get a nicer behavior than in 
Fig | 
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A variable xmin but Q=2 
S(pt) with pto=0.92 
A variable xmin and Q 
S(pt) with pto=0.64 
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Figure 17: The screening function S(pj_) fitted to A. Our parameterization of the screening 
function can best be fitted to the curve with a fix Q. 



3 
2.5 
2 



1 

0.5 




I I I II, 

1. 




' 1 


2. 






3. 

088(1 x,..,.)° 075 
- 0.59(1 x,..,.) 005 " 
0.4(1/x min ) 058 


■ 


■ 


■ 






■ 






- A 

, I , , I 


, I 


. 1 



100 



1000 



1/x 



10000 



100000 



min 



Figure 18: How p±o varies with the energy for different start configurations. Scenarios 
1,2 and 3 according to the notation in the text. 
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4.5 The Radius of the Proton in our Model 



As another application of our model we are going to look at how the radius of our model 
of the proton evolves as we increase the energy of the system. 
We do this by plotting the parton density function f(r): 



f(r) 



1 dn. 



partons 



dn, 



dr 



partons 

d 2 r 



(37) 



the phase-space factor 1/r is introduced to avoid having a vanishing distribution at r = 0. 
The center, from which the distance r is calculated, is simply the total average of the 
position of all partons. 
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Figure 19: How the radius of our proton model varies with the energy. 



By defining some lowest density of the proton it is easy to define a radius as the 
distance from origo to the point on the curve at this lowest density. To vary the energy 
we once again changes the parameter x min , see Fig |19[ Here the density function is plotted 
for some different x min . As is expected, the radius increases with the energy. We see that 



the curves in Fig 19 are almost parallel to each other. This means that the partons stays 



rather close to their original parton mother. Though a closer examination of Fig [19] shows 
that the curves spreads towards larger R. This is an indication of some diffusion, i.e. 
partons that have traveled far from their original position. 
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5 Summary and Conclusions 



This paper is divided in two main parts: 

1. The Proton Model 

2. The Studies of the Model 

In the first part we describe how the model is constructed and which conditions we have 
made upon it. 

The model is iterative in its structure; we start with a configuration of only 5 partons. 
These are then evolved towards higher virtuality scales Q to some Q ma x where we stop 
the evolution. The evolution is done in an exclusive picture i.e. we follow each parton by 
simulating the probability for a branching by random numbers. To know the probability 
for a branching we use the Altarelli-Parisi splitting kernels. 

If we allowed all partons there would be infinitely many in our model, and this is not 
possible. To avoid this we introduced a cutoff in momentum space. The branchings that 
produces partons with a momentum below this x m i n were not allowed, all the partons are 
still kept into the parton shower though. Despite this rather abrupt cutoff, the distribution 
of the partons is not changed significantly. 

We have also made the model so that it is always colour neutral. This is done by 
having a colour neutral start configuration. Then, as we evolve, the rules of QCD are 
followed and no partons are thrown away. 

All partons were assigned a transverse position in the coordinate space and a longitu- 
dinal value of the momentum. This is done to try to avoid difficulties from the Heisenberg 
uncertainty principle. 

We have compared our model to the model by the GRV group, see Fig. || to Fig. 0. 
Our model cannot compete with GRV but we have a rather good agreement between the 
two models, especially at high Q. 

In the second part of this paper we studied the so called screening effect with our 
model. The screening effect is that the partons will screen each other so the charge seen 
from outside the proton will be lower than the naive incoherent sum. By using our model 
as a target we sent gluons with different momentum at the proton. 

We made some different assumptions of how the gluon interact with the proton model: 

• A running Q: That the incoming gluon will react mostly with partons with Q (and 
thereby p±) equal to the gluon momentum 

• A running x m i n : The cutoff variable x m i n should vary with the Q scale. 

• The distance between the partons: When a branching occurs the distance between 
the daughters are decided by Eq (p4|) , this is varied 

The different curves of this is presented in Fig. |18| We see that there are some 
differences, but the qualitative behavior is the same. We clearly have the screening effect 
and the a jet will be better behaved than in Fig. 
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